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Abstract 

A prototype Silicon- Tungsten electromagnetic calorimeter (ECAL) for an International 
Linear Collider (ILC) detector was installed and tested during summer and autumn 
2006 at CERN. The detector had 6480 silicon pads of dimension 1x1 cm^. Data were 
collected with electron beams in the energy range 6 to 45 GeV. The analysis described 
in this paper focuses on electromagnetic shower reconstruction and characterises the 
ECAL response to electrons in terms of energy resolution and linearity. The detector 
is linear to within approximately the 1% level and has a relative energy resolution of 
(16.6 ± 0.1)/-/B(GeV) e 1.1 ± 0.1 (%). The spatial uniformity and the time stability of 
the ECAL are also addressed. 

Key words: 

CALICE, ILC, electromagnetic calorimeter, silicon detector, electron reconstruction 



Preprint submitted to NIM A 



November I4, 2008 



1. Introduction 



The CALICE Collaboration is conducting R&D into calorimctric systems for the 
ILC [H — a proposed e+e" linear collider intended to operate at a centre of mass energy 
ranging up to the TeV scale. The physics scope at the ILC includes precise measurements 
of the triple- and quartic-gauge bosons interactions, as well as the complete character- 
isation of the Higgs and top quark sectors. In addition, hints of physics beyond the 
Standard Model could be addressed in a model-independent way. 

The final states are typically multiple hadronic jets, accompanied frequently by low 
momentum leptons and/or missing energy. In such cases, lepton identification is difficult 
and the signature of the final states of interest relies on the identification of Z or/and W 
bosons in their decay modes into two jets. In order to distinguish them efficiently, a jet 
energy resolution close to 30%/-\/i?/GeV has to be achieved A precise reconstruction 
of the jet direction is also required. These are the main requirements driving the detector 
design in general at the ILC and the calorimetry design in particular. 

The target jet energy resolution represents an improvement by a factor of two over the 
best obtained in previous detectors. Moreover the detection environment becomes more 
complex with increasing centre of mass energy. A promising way to achieve this increase 
in resolution is through designing a detector system optimised for the so called "particle 
flow" approach which relies on the separate reconstruction of as many particles in 
the jet as possible, using the most suitable detector systems. 

The success of such an algorithm depends on the quality of the pattern recognition in 
the calorimeters. For particle flow, a high spatial granularity is therefore as important as 
the intrinsic energy resolution for single particles. Furthermore, the overall design of the 
detector (tracking, electromagnetic and hadronic calorimetry) needs to be considered in 
a coherent way. 

The design of the ILC detectors can be optimised using Monte Carlo simulations, but 
in order to do this, it is crucial to validate the Monte Carlo tools with data. Therefore, 
the R&D of the CALICE Collaboration has two broad aims. The first is to construct 
realistic calorimeter prototypes, and learn about their operation and behaviour in beam 
tests. The second objective is to compare the data with Monte Carlo simulations using 
the same tools used for the full detector. This is especially important in the case of 
hadronic showers, where many models are available, which make differing predictions for 
the calorimeter response. The CALICE plan is to expose complete calorimeter systems 
(electromagnetic and hadronic, using various technologies) to test beams of electrons, 
muons and hadrons. To this end, a first round of beam tests was performed at DESY 
and CERN in summer 2006, using a Silicon- Tungsten sampling electromagnetic calorime- 
ter followed by a hadron calorimeter composed of iron and scintillator tiles [1], and 
then a Tail Catcher and Muon Counter (TCMT) of iron instrumented with scintillator 
strips 0]. 

In this paper, we report results of exposure of the prototype to electron beams in 
the energy range 6-45 GeV at the CERN 116 beam line . In Section [2] we outline 
the layout of the beam tests. The ECAL is briefly described in Section |3] and some key 
technical aspects of its performance are highlighted. Section |4] summarises the Monte 
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Carlo simulation. Features of the electron beam data are reviewed in Section [5] and 
the uniformity across the detector is addressed. The results of the energy measurement 
together with some of their systematic uncertainties are presented in Section [S] for the 
detector areas of uniform response. 

Several other studies of the prototype are ongoing, exploiting its unprecedentedly fine 
segmentation and capacity to observe shower development in detail. These studies will 
be reported in subsequent publications. 



2. Experimental Setup 

A sketch of the CERN H6 Q test beam setup is presented in Figure [T] and a detailed 
description of the detectors can be found in [3|. The coordinate system used is right 
handed. The surface of the drift chamber (DCl) closest to the ECAL defines the origin, 
the z axis is the beam axis and x and y the horizontal and upward- vertical, respectively. 

The physics program and the overall electron, pion and muon statistics collected are 
extensively discussed in The beam trigger was defined by the coincidence signal of 
two scintillator counters. In addition, three drift chambers were used to monitor the 
beam. A threshold Cerenkov detector was also available for c/tt discrimination. 

This paper presents the measurement of the ECAL response to electrons normally 
incident on the calorimeter surface. The event display for one of these events is shown 
in Figure [21 where the energy of the hits is measured in MIP units, one MIP being the 
signal left by a minimum ionising particle. 
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Figure 1: Sketch of the CERN test beam setup. The the right handed coordinate system used hereafter 
is indicated. 



3. The Si-W ECAL prototype 

A detailed description of the ECAL hardware is given in [3], along with details of 
the commissioning and a number of technical features of the system calibration and 
performance. The ECAL prototype consisted of 30 layers of tungsten, the first ten 
of thickness 1.4 mm, the next ten of 2.8 mm and the last ten of 4.2 mm, comprising 
24 radiation lengths in total at normal incidence. The development of the showers was 
sampled using 30 layers of silicon PIN diode pads interleaved between the tungsten plates. 
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The silicon thickness was 525 /im, with each pad having a transverse area of 1x1 cm^. 
The sensors were implemented on 4 square inch wafers, segmented into 6x6 pads. At 
the time of the 2006 CERN beam tests, each layer consisted of a 3x2 array of wafers, 
i.e. 18 pads horizontally and 12 pads vertically, leading to a total of 6480 pads for the 
ECAL. 

Blocks of random triggers were recorded during data taking in order to monitor 
pedestals and noise. Short term changes and shifts in pedestals caused by large sig- 
nals in neighbouring cells were monitored and corrected using cells without signal in 
beam events [3]. The uncertainty on the pedestal levels was estimated to be less than 
0.002 MIPs, negligible compared to the energy deposited by electron showers (a 10 GeV 
electron deposits on average 1450 MIPs). The noise level was typically 0.13 MIPs; its 
spread channel-to-channel was 9% of the mean noise and the spread run-to-run was less 
than 1% of the mean noise. The low spread of the noise justifies the use of a single energy 
threshold for all cells in the detector. 

Calibration constants for each pad were determined using muon events. The response 
of each cell was fitted by a convolution of a Landau distribution with a Gaussian. The 
most probable value of the underlying Landau function was taken to define the MIP value 
for each cell, and the raw energy for each cell in data was corrected to units of MIPs. 
All but 9 pads were functioning and successfully calibrated. The calibration constants 
were determined with an accuracy of 0.5% and had a cell to cell dispersion of 5%. Data 
taken in the various beam test periods during summer and autumn 2006 showed well 
correlated calibration constants, with differences less than 1.6%. 

One feature of the data which has not been accounted for in the detector simulation is 
associated with showers which deposit a sizeable energy in the guard ring surrounding a 
wafer. This is a cause of correlated crosstalk, observed as a distinctive square pattern of 
low energy hits in a number of cells around the periphery of the wafer. The prevalence of 
this effect increases with the shower energy crossing the guard ring and therefore its rate 
is significantly reduced when considering only electrons impinging on the centre of the 
wafers. In the future, the design of the guard rings will be modified in order to prevent 
this problem. 

After calibration, the ECAL data consist of hits in the cells of the calorimeter with 
energies in units of MIPs. In order to remove most of the noise signals, a threshold cut 
of 0.6 MIP was imposed on each cell, almost five times the mean noise level. 

4. Monte Carlo Simulation 

The test beam setup is simulated with Mokka a Geant4 0]-based Monte Carlo 
program, followed by a digitisation module simulating the response of the data acquisition 
electronics. The material upstream of the ECAL is taken into account. The subdetectors 
are simulated with different levels of detail, depending on their impact on the physics 
analysis: material simulation only for the Cerenkov detectors, raw energy depositions 
stored for the trigger counters, partial electronics simulation for the tracking detectors. In 
the case of the ECAL, the simulation gives the raw energy depositions in the Si pads and 
the readout electronics is simulated assuming that each channel exhibits only Gaussian 
noise. The mean values of the noise for each channel follow a Gaussian distribution over 
the ECAL, with the mean value of 0.13 MIPs and 0.012 MIPs dispersion, as measured. 
Since 99.9% of the ECAL cells were functioning, the impact of the non-responding cells 
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is expected to be small and their signals were not supressed in the simulation for this 
analysis. 

The beam simulation assumes a parallel beam with Gaussian width reproducing the 
observed beam profile. To study systematic effects due to lateral leakage of the showers, 
samples are also generated with a beam spread uniformily over the ECAL front face. A 
Gaussian momentum dispersion consistent with the settings of the beam collimators Q 
is applied for each run. 



5. Selection of Electron Events 

Single electron showers are selected using the energy recorded in the ECAL. This 
energy, £^raw, is calculated with the three ECAL modules weighted in proportion to the 
tungsten thickness: 

i=9 i=19 i=29 

^^raw = ^ + 2 ^ + 3 ^ £;„ (i) 

i=0 1=10 4=20 

where Ei is the energy deposit in layer i. The distribution of i?raw is shown in Figure [3] 
for a typical 15 GeV event sample. The electron peak at around 3900 MIPs is clearly 
visible; however, the muon and pion contamination in the beam gives an additional peak 
at 85 MIPs and the region between the two main peaks is populated with pions. Electron 
candidates are selected by requiring: 

, graw(MIP) 
< i?.ea.(GeV) < 

The significant pion contamination present in some of the data runs is reduced by de- 
manding a trigger signal from the threshold Cerenkov counter in the beam. The effect 
of this additional requirement is indicated by the shaded region in Figure [3] 

5.1. Rejection of the beam halo 

The rejection of the beam halo is implemented run-by-run. The x and y acceptance 
for the incoming electron track is chosen such as to achieve a reasonably flat distribution 
of the mean energy deposition in the ECAL. 

5.2. Inter-wafer gap effect 

Around the pads in each wafer, a non-active region of 1 mm width was used for 
a grounded guard ring structure. This creates a non-active gap between adjacent Si 
pads situated on different wafers (2 mm) which is significant compared to the transverse 
shower size. These non-active regions, called in the following "inter-wafer gaps" degrade 
the prototype response when showers traverse them. This is illustrated in Figure |4] for 
30 GeV electrons impinging on the calorimeter at normal incidence. Here the mean value 
of -Eraw (Equation [1]) is plotted as a function of the shower barycentre {x,y), defined as : 

=^(E;,x„£;,2/,)/I]^^4 (3) 



The sums run over all hit cells in the calorimeter. Dips in response corresponding to 
the guard ring positions are clearly visible: the energy loss is about 15% when electrons 
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Figure 3: Distribution of total ECAL hit energies for a 15 GeV electron run with a significant pion 
content. The -Braw selection window and the shaded area obtained by demanding a signal from the 
Cerenkov counter are shown. 



impinge in the centre of the x gaps and about 20% in the case of the y gap. In order 
to recover this loss and to have a more uniform calorimeter response, a simple method 
was investigated. The ECAL energy response, f{x,y) = i^raw/^'boam, is measured using 
a combined sample of 10, 15 and 20 GeV electrons, equally populated. 

The response function /, normalised such as to have a unit response in the middle of 
the wafers, is displayed on Figure [3 To characterise the x (y) response, the events were 
required to be outside the inter- wafer gap in y (x) , leading to an important difference in 
the number of events for the two distributions, since the beam is centred on the y gap. 
It can be parametrised with Gaussian functions, independently in x and y: 

m. p) = (l - ..exp (-^i^)) (l - «..xp (- ) (4) 

Here, Xgap and ?/gap are the positions at the centres of the gaps in x and y, respectively, 
ttx (ciy) and ax (ay) their respective depths and widths in the two directions. The results 
of the Gaussian parametrisations are given in Table [TJ Thegap in x is shallower and 
wider than that in y, due to the staggering of the gaps in a; [3[. 





position (mm) 


(T (mm) 


a 


X direction 


-30.0 


4.3 


0.143 


y direction 


-8.4 


3.2 


0.198 



Table 1: Gaussian paramctrisation of the inter- wafer gaps. 



As illustrated in Figure [SI when the energy of each shower is corrected by 1//, the 
average energy loss in the gaps is reduced to a few percent level. The low energy tail in 
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Figure 4: Mean values of Eraw for 15 GeV electrons as a function shower barycentre, transverse to the 
beam direction. The energies have been scaled down by a factor 266 to provide approximate conversion 
to GeV. 
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Figure 5: Normalised f{x,y) as a function of the shower barycentre coordinates, for a combined sample 
of 10, 15 and 20 GeV electrons. 
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the energy distribution is also much reduced (Figure [7]). The correction method rehes 
only on calorimetric information and can be applied both for photons and electrons. 
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Figure 6: Mean -Braw as s function of the shower barycentre coordinates for 20 GeV electrons, before 
(open triangles) and after the inter-wafer gap corrections (solid circles) were applied on -C/raw ■ 



Even though it is possible to correct for the inter-wafer gaps on average, for individual 
events their presence will induce fluctuations in the energy response and degrade the 
ECAL resolution compared to a continuous calorimeter. In the data described here, the 
beam centre was close to the inter-wafer gap in y, artiflcially increasing the impact of 
the inter-wafer gaps compared to an experiment with a beam uniformely spread over the 
ECAL front face. Moreover, since the beam width varies strongly with the beam energy, 
the impact of the energy lost in the gaps is different at each energy. Therefore, in order 
to assess the energy response and resolution of the prototype as a function of energy in a 
unbiased way, only particles impinging in the middle of the wafers are selected. Since the 
gap effect is negligible more than 4 standard deviations away from the gap centre, the 
shower barycentre for selected events is required to be at a distance larger than 17.2 mm 
from the centre of the inter- wafer gap along x and 12.76 mm away from the centre of the 
y gap. 

A sustained R&D effort is being made to reduce the non-active areas, both by reducing 
the size of the inter-wafer gaps and by increasing the size of the wafers. The next Si- 
W ECAL prototype will have 9 x 9 Si pads in a wafer which leads to a significant decrease 
of the non-active areas. 



5.3. Selection of showers well contained in the ECAL 

The fiducial volume in which the showers are fully contained in the ECAL was esti- 
mated using electrons away from the inter-wafer gaps and pointing at the centre of the 
ECAL. The radial shape of an average 45 GeV electron shower is shown in Figure [3 
both for data and simulation. The simulation reproduces the shower width to better 
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Figure 7: Energy distribution for 20 GeV electrons in the cases of: events outside the inter-wafer gaps 
(sohd histogram), all events without inter- wafer gap corrections (open histogram) and all events with 
inter- wafer gap corrections (solid circles). The histograms are normalised to the same number of entries. 

than 2%: 95% of the shower energy is contained within 30.5 mm {i.e. less than four 
pads), to be compared with 29.9 mm in the case of the simulated showers. To ensure 
radial containment, all electrons impinging on the ECAL front face less than 32 mm from 
one of the ECAL borders are therefore excluded from the selected sample. 

The longitudinal containment of the showers is ensured by rejecting events which 
have the maximum of the energy deposited along the z direction in the first five layers 
or the last five layers of the prototype. Only 0.21% of the simulated 6 GeV electrons fail 
these containment criteria and 0.02% of the 45 GeV electrons. 

5.4- Rejection of electrons showering in front of ECAL 

The data recorded at CERN contain a significant number of events which have ap- 
proximately the expected energy for a single electron, but whose spatial structure clearly 
exhibits double showers. A likely explanation is bremsstrahlung far upstream in the 
beamline. In the Monte Carlo simulation, the known material between the Cerenkov 
counter and the calorimeter is simulated, and yet the agreement between the rate of 
double shower events is poor between simulations and data. Before comparing data and 
Monte Carlo, it is therefore necessary to select a sample of single electron showers. To 
this end, the energy deposits in the shower are projected in a two-dimensional histogram, 
on the transverse, x — y, plane. The binning of the histogram is the same in x and y and 
corresponds to the cell size (1 cm). A simple nearest- neighbour clustering algorithm (in- 
cluding diagonal neighbours) is applied on the bins with energies above a given threshold 
T, in order to select events with more than one local maximum for the energy deposit. 
For each event we determine the maximum value of the threshold, Tmax, above which the 
event would be reconstructed as a single cluster. In Figure Owe compare the distribu- 
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Figure 8: Energy deposited in ECAL as a function of the radial distance to the longitudinal shower axis, 
integrated over 171671 showers of 45 GeV. 



tions of Tinax between data and simulation for a 30 GeV beam. A sizeable discrepancy is 
seen for larger values of Tmax and therefore a cut is applied on Tmax- The cut is energy 
dependent, varying from 50 MIPs at 10 GeV to 120 MIPs at 45 GeV. This cut typically 
rejects ^20% of data and '^2-3% of simulated events at the higher energies. 

A summary of the selected electron and positron data is shown in Table [21 The 
number of simulated events available for each energy is also indicated. 



Energy (GeV) 


date 


data statistics (kevts) 


MC statistics (kevts) 


6 


Oct 


6.6 


83.2 


10 


Aug, Oct 


43.1 


80.3 


12 


Oct 


27.2 


72.8 


15 


Aug, Oct 


51.4 


70.3 


20 


Aug 


62.9 


56.2 


30 


Aug 


42.3 


55.2 


40 


Aug 


22.9 


67.8 


45 


Aug 


108.6 


108.8 



Table 2: Summary of the electron events selected for this analysis. 
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Figure 9: Distribution of the variable Tmax (described in the text), which is used to reduce the contri- 
bution of double showers. Data and simulation are compared at 30 GeV. 

6. Performance Studies 

6.1. ECAL Sampling Fraction Scheme 

The ECAL is made of 30 layers grouped in three modules of 10 layers each Q. Each 
tungsten sheet has the same thickness in a given module. However, as can be derived 
from Figure I10[ where one passive tungsten layer sandwiched between two active silicon 
layers is shown, two successive silicon layers are either separated by one thickness of 
tungsten or by the same thickness of tungsten plus two thicknesses of PCB, aluminium 
and carbon-fibre-epoxy composite. A different sampling fraction, defined as the ratio of 
the energy deposited in the active medium to the total energy deposit (sum of the energy 
deposits in the active and passive medium), is therefore expected for the even and the 
odd layers of the same calorimeter module. 

The easiest method to investigate this difference is to compare in each module the 
mean energy deposits in odd and even layers. For the first module, if we neglect the 
shower profile, the ratio of the two is 



with rj being, approximately, the ratio of the non-tungsten radiation length to the tung- 
sten radiation length. 

When counting the layers starting from zero, the odd layers are systematically shifted 
compared to the even layers towards the shower maximum and the measurement of R is 
biased by the shower development. To overcome this bias, R is measured twice, either 
comparing the odd layers with the average of the surrounding even layers, or comparing 



R = 



^evcn 



1 + ^, 



(5) 
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Figure 10: Details of one ECAL slab, showing one passive tungsten layer sandwiched between two active 
silicon layers. The dimensions are in mm. In contrast to the upper silicon layer which is preceded by a 
layer of tungsten only, the lower silicon layer is preceded by a larger passive layer: the PCB, aluminium, 
glue and carbon structure as well as the tungsten. 



the even layers with the average of the neighbouring odd layers: 

^, ^ (El + E3 + E5 + E7) 

/ Ea+E2 I E2+E4 I E4+E6 I Eg+Es \ 
\2^2'2'2/ 

/ E1+E3 I E3+E5 I E5+E7 , Et+E9 \ 

j^ll _ \ 2~'~2'2'2/ 

(E2 + E4 + Eg + Eg) 

where E„ is the energy deposit in the layer number n and the brackets indicate that 
mean values are used. The value of 77 is taken as the average of i?' — 1 and R" — 1, while 
the difference between them gives a conservative estimate of the systematic uncertainty 
due to the shower shape. As an example, the distributions of the energy deposits in the 
odd and even layers are shown in Figure fTTl for 20 GeV electrons. 

The values of rj, obtained using the first module and for different beam energies, are 
displayed in Figure [121 The overall value is (7.2 ± 0.2 ± 1.7)%. The measurement of rj 
using the second and third module gives compatible results and the corresponding value 
obtained from simulation is (4.7 ± 0.2 ± 2.0)%. 

In computing the total response of the calorimeter, the sampling fraction for layer i 
is given hy Wi ~ K for even layers and Wi = K + rj for the odd layers, with K = 1, 2, 3 
in modules 1, 2, 3, respectively. 

6.2. Linearity and energy resolution 

The total response of the calorimeter is calculated as 

£;,ec(MIPs) = ^u;,£;, (8) 

i 

with Wi the sampling fraction for the layer i. Its distribution for electrons at 30 GeV is 
shown in Figure [T51 together with a fit using a Gaussian function in the range [—a, +2a]. 
There is reasonably good agreement between data and simulation. An asymmetric range 
is chosen for the fit in order to reduce sensitivity to pion background, to radiative effects 
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(6) 
(7) 




Energy (MIPs) 

Figure 11: Energy deposits in odd layers (continuous histogram) and average energy deposits in their 
surrounding even layers (dashes) by 20 GeV electrons, in the first ECAL module. 
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Figure 12: Values of 77 as a function of the beam energy. The uncertainties are statistical and the dashed 
line gives the average value of rj. 
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upstream of the calorimeter, and to any residual influence of the inter-wafer gaps. The 
position of the peak is the mean energy response (called in the following -Emoan) and its 
distribution is shown in Figure [TJ] as a function of the beam energy. The uncertainties 
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Figure 13: Gaussian parametrisation of -Broc for 30 GeV electron data (dashes). The range of the fit is 
[— cr, +2(t]. The data (solid squares) and simulation (continuous line) arc superposed. 



From the dispersion of -Emcan in the different runs at the same nominal beam energy, 
the uncertainty of the beam mean energy, -Bboam, was estimated to be 



0.12 



Ehc 



Sbcam(GeV) 



®0.1%, 



(9) 



The first term is related to hysteresis in the bending magnets, while the calibration and 
the uncertainties on the collimator geometry give the constant term. For comparison, 
in the uncertainty on the beam mean energy was quoted as 



0.25 



Ehc 



Sbcam(GeV) 



® 0.5%, 



(10) 



The first of these parametrisations of the uncertainty (Equationl^l) is used in the following, 
except for checks of systematic uncertainties. 

The mean energy response can be parametrised as -Bmcan = P • ^-bcam — ck, while the 
measured energy iJmoas is given by i^mcas = ^-mcan + ct- The parameter is a global MIP 
to GeV calibration factor. The offset a is partly due to the rejection of the low energy 
hits and it increases steadily with the hit energy threshold, as displayed on Figure [151 
On the same figure are also shown the values of the offset, as expected from simulation. 
The bias introduced by the fact that the uncertainty on the beam mean energy decreases 
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Figure 14: Energy response of the ECAL as a function of the beam energy. For clarity, all the runs 
around the same nominal energy of the beam were combined in one entry for the plot, for which entry 
the uncertainty was estimated assuming that the uncertainties on the individual runs were uncorrelated. 



with increasing beam energy and enhances therefore the weight of the high energy runs 
was estimated by artificially assigning to the simulated data uncertainties on i?boam 
according to Equation [51 By taking into account this bias, the disagreement between 
data and simulation is somewhat reduced (Figure [TS]). 

The residuals to linearity of the measured energy, converted from MIPs to GeV using 
a constant 266 MIPs / GeV conversion factor (obtained from Figure [T4|l , are shown in 
Figure [12] as a function of the beam energy. The residuals are within approximately the 
1% level and are consistent with zero non-linearity. Data and simulation agree within 
one standard deviation. 

The relative energy resolution, Ai^mcas/^'mcas , as shown in FigurefTTl can be parametrised 
by a quadrature sum of stochastic and constant terms 

where the intrinsic momentum spread of the beam was subtracted from the ECAL 
data By varying the range for fitting i?rec, the variation of the stochastic term 
has the same order of magnitude as the statistical error, whereas the constant term 
remains stable: a fitting interval reduced to [— 0.75ct, +1.75ct] improves the stochastic 
term to 16.5 ± 0.2, whereas an enlargement to [— 2ct, +3cr] degrades it to 16.7 ± 0.1. The 
expected resolution from simulation is 

/l7£±0T ^^^^^ \ 

B(GeV) 7 
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Figure 15: Variation of the linearity offset with the hit energy threshold: data (solid squares), data with 
the fitting range of Erec enlarged to [—2a, +3cr] and considering the parametrisation in Equation [10] for 
the uncertainty on the beam mean energy (solid triangles), Monte Carlo (open squares), Monte Carlo 
with an artificial uncertainty on the beam mean energy as given by Equation 1101 (solid circles) . For 
clarity, the points were artificially shifted along the x axis around the nominal -EhH threshold. 
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Figure 16: Residuals to linearity of Smeas as a function of the beam energy, for data (solid squares) and 
simulation (open squares). All the runs around the same nominal energy of the beam were combined 
in one entry, for which the uncertainty was estimated assuming that the uncertainties on the individual 
runs were uncorrelated. For clarity, the Monte Carlo points were artificially shifted along the x axis 
around the nominal -Ebcam- 
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which agrees within 5% with the measured resolution of the prototype. 
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Figure 17: Relative energy resolution (Ai?niGas/i?mcas) as a function of the beam energy (solid squares), 
and its usual parametrisation as s/ \/E © c. For clarity, the 35 runs available were combined into 8 
different beam energy points for the plot. For the parametrisation of the energy resolution each run was 
however treated individually. The values expected from simulation are superposed (open squares). 



Different systematic checks have been performed on the data. Variations of the hnear- 
ity and resolution against the minimal accepted distance between the shower barycentre 
and the nearest inter-wafer gap, when the energy threshold for considering the hits is 
0.6 MIPs are shown in Table [H In addition, this hit energy threshold has itself been 
varied (Table[3]). In order to investigate the potential effects linked to the beam position, 
the energy response is also compared for showers with barycentres located in the right 
hand side (negative x coordinates) and in the upper half of the detector (upper row of 
wafers) as summarised in Table [S] The results of all checks are consistent. Since data 
were taken in both August and October 2006, it was also possible to check the response 
stability in time and no significant differences between the two data samples are observed. 



7. Conclusion 

The response to normally incident electrons of the CALICE Si-W electromagnetic 
calorimeter was measured for energies between 6 and 45 GeV, using the data recorded 
in 2006 at CERN. 

The calorimeter response is linear to within approximately 1%. The energy resolution 
has a stochastic term of (16.6 ± 0.1)%/y^?(GeV), whereas the constant term is 1.1 ± 
0.1%. Several sources of systematic uncertainties have been investigated and their effect 
is within the statistical uncertainties. The agreement between data and Monte Carlo 
simulation is within 5%. 
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shower distance to the gaps (in standard deviations) 




A 






(linearity) 


16.8/32 


17.6/32 


18.9/32 


24.2/32 


Oi 

(MIPs) 


93.9±11.1 


96.3±11.2 


97.8±11.5 


99.1±11.6 


(MIPs/GeV) 


266.3±0.5 


266.6±0.5 


266.8±0.5 


266.8±0.5 


resolution (%) 
(stochastic term) 


16.7±0.1 


16.6±0.1 


16.4±0.2 


16.3±0.2 


resolution (%) 
(constant term) 


l.OiO.l 


l.OibO.l 


l.liO.l 


1.2±0.1 



Table 3: Impact of the distance of shower to the intcr-wafer gaps on the ECAL Unearity and resolution. 
The distance is given in terms of standard deviations to the gap centre, with the standard deviation 
defined by the Gaussian pajrametrisation of the gaps. 





E^^■^^ cutoff (MIPs) 




0.5 


0.7 


0.9 


X^lndf 
(linearity) 


18.0/32 


17.8/32 


18.0/32 


a 


93.0 ±11.2 


98.9±11.1 


105.6±11.1 


(MIPs) 








/3 


266.8±0.5 


266.3±0.5 


265.8±0.5 


(MIPs/GeV) 








resolution (%) 


16.6±0.1 


16.5±0.1 


16.6±0.1 


(stochastic term) 








resolution (%) 


l.OiO.l 


l.libO.l 


l.liO.l 


(constant term) 









Table 4: Impact of the hit energy cutoff on the ECAL linearity and resolution. 
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right side 


upper part 


a (MIPs) 


96.1±10.9 


97.7±11 


(3 (MIPs/GeV) 


266.6±0.5 


266.8±0.5 


resolution (stochastic term) (%) 


16.8±0.1 


16.8±0.2 


resolution (constant term) (%) 


l.liO.l 


l.liO.l 



Table 5; Response to electrons crossing the right hand side and the upper part of the ECAL. 
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